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ABSTRACT 

Observational evidence for the radial alignment of satellites with their dark matter host has been 
accumulating steadily in the past few years. The effect is seen over a wide range of scales, from 
massive clusters of galaxies down to galaxy-sized systems, yet the underlying physical mechanism has 
still not been established. To this end, we have carried out a detailed analysis of the shapes and 
orientations of dark matter substructures in high-resolution N-body cosmological simulations. We 
find a strong tendency for radial alignment of the substructure with its host halo: the distribution 
of halo major axes is very anisotropic, with the majority pointing towards the center of mass of the 
host. The alignment peaks once the sub-halo has passed the virial radius of the host for the first time, 
but is not subsequently diluted, even after the halos have gone through as many as four pericentric 
passages. This evidence points to the existence of a very rapid dynamical mechanism acting on these 
systems and we argue that tidal torquing throughout their orbits is the most likely candidate. 
Subject headings: methods; N-body simulations, galaxies: kinematics and dynamics, galaxies: clusters 
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1. INTRODUCTION 

Anisotropy in galaxy orientations has been a matter 
of debate for several decades, and many conflicting re- 
ports have been published. Past studies have found ev- 
idence for three dif ferent types of al ignment: alignment 
betw een clusters (jBinggelil Il982l : iPlionis fc Basilakoi 
I2002D . between the brightest clust er galaxy ( BCG) 
and the satellite distribution (e.g., I Yang et al] 120061 ) 
and between the orientation of satellites and their 
host (iHawlev fc Peebleal Wnk iPereira fc KiThl] [20051 : 
I Agustsson fc BrainerdI I2006D . This last type of align- 
ment, which we will refer to as radial alignment and is 
the focus of this paper, has been the h ardest to confirm 
(jTrevese et al.lll992l : [Torlina et al.ll2007[ ). since it requires 
high quality data on small scales. In recent years, this 
field has seen a resurgence, largely due to the arrival of 
the S loan Digital Sky Survey (SDSS) (jAbazaiian et al.1 
l2005f ). The SDSS provides accurate measurements of 
isophotal shapes for millions of galaxies, and this has 
finally allowed large st atistical studies of gala xy align- 
ments to be performed. IPereira fc KuhnI (|2005l ) targeted 
galaxies in massive X-ray selected clusters and found 
a significant tendency for their radial alignment . This 
result has since been confirmed by Faltenbacher et al.l 
(|2007( ) for a larger sample of gr oups optically selected 
from t he SDSS. On smaller scales. IXgustsson fc BraineTdl 
()2006f ) found a tendency for satellite galaxies in the SDSS 
to be radially aligned with their host g alaxy, whereas on 
large scales. iMandelbaum et al.l (|2006a ) found a very sig- 
nificant correlation between the orientations of galaxies 
and the surrounding density field traced by galaxy over- 
densities. 

Initially motivated by the prospect of using galaxy ori- 
entations to probe their formation histories, these studies 
are now also driven by the need to calibrate weak lens- 
ing and cosmic shear measurements. A key assumption 
for lensing techniques is that the population of galaxies 
being lensed is randomly oriented. Some intrinsic align- 
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ments between galaxies can be dealt with easily: e.g. 
downweighting close pairs readily removes contamination 
by alignments indu ced in interacting systems. However, 
as lHirata fc SeljaS (2003) pointed out, if galaxy orienta- 
tions are affected by their surrounding density field {e.g. 
a galaxy cluster), then they will also be correlated with 
the orientations of the background population of galax- 
ies that is being lensed by that field. This correlation 
between widely separated redshift bins cannot trivially 
be removed. 

Given the growing body of evidence suggesting that 
galaxy orientations are anisotropic, and the pressing need 
for an accurate quantification of intrinsic alignments for 
weak lensing, it seems crucial that we try and find the 
physical cause behind these anisotropics. There are two 
commonly pr oposed explanations. The first, initially 
developed by iPeebled |l969) in his tidal torque theory 
(TTT), explains the anisotropy as a left-over primor- 
dial effect. TTT ascribes the orientation and rotation 
of galaxies to torquing during their formation. It there- 
fore follows that the signal should be stronger on the 
outskirts of the cluster, and that it wanes with time, 
such that older, more relaxed clusters should exhibit less 
te ndency for alignm e nt. T he other alternative, proposed 
bv IPereira fc KuhnI (120051 ). is a dynamical mechanism, 
i.e. an interaction with the tidal field of the host clus- 
ter that gets progressively stronger during infall and is 
not erased by subsequent orbital motions. Observational 
studies have so far been unable to distinguish between 
the two, due to difficulties in constraining galaxy orbits 
and in measuring galaxy shapes accurately out to large 
redshifts. 

A different approach is needed, and a few numeri- 
cal studies have recently been published on this sub- 
ject. Studies of simulated halo shapes and orientations 
have been performed around voids (Brunino et al. 2007), 



along filanients and sheets dAragon-Calvo et al.l 1200 
lAltav et al.ll2006l: iHahn et al.l l2007D and in a Milky- Way 
tvpehalo (jKuhlen et al.l2007| ). Anisotropies are reported 
in all three environments. The advantages of working 
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with simulated clusters are obvious - 3D spatial infor- 
mation means we do not suffer dilution from projection 
effects. Also, with enough temporal and spatial reso- 
lution, we can follow the galaxies as they fall into the 
cluster along filaments, and beyond, as they orbit inside 
the cluster. By tracking the effect's evolution with time, 
we will be able to more precisely determine its source. 

We start (§2.1) by introducing the simulations used for 
this analysis and describing the properties of the eight 
host halos. Our methods for finding the substructure 
halos (§2.2) and determining their shapes (§2.3) follow, 
along with a study of the reliability of our shape measure- 
ments. With this information, we then show in section 
3.1 that cosmological dark-matter simulations do indeed 
produce radial alignment in clusters at z w 0, and we 
study the correlation of this effect with various parame- 
ters, such as host halo mass and distance from the clus- 
ter center (§3.2). Having established the importance of 
the alignment effect in our simulations, we use the high 
temporal resolution to study its evolution with time in 
section 3.3, and its dependence on orbital phase (§3.4). 
We argue in section 4.1 that tidal torquing by the host 
halo tidal field is responsible for the alignment of sub- 
structure, and compare our results with previous obser- 
vations in §4.2. We end (§4.3) by briefly speculating on 
the possible consequences of such a mechanism for the 
morphological and orbital evolution of galaxies in clus- 
ters. 

2. SIMULATIONS AND ANALYSIS 

2.1. The Data 

The iV-body sim ulations us e d in this work are pre- 
sented in detail in I Gill et ahl (j2004[ ). and we describe 
them here only briefly. Using t he open source a daptive 
mesh refinement code MLAPM (K nebe et al.ll200lh . a set 
of four initial conditions at redshift z = 45 in a standard 
ACDM cosmology (r^o = 0.3, ilx = 0.7, ilbh^ = 0.04, h = 
0.7, (78 = 0.9) were created. From an initial distribution 
of 512^ particles in a box 64:h~^ Mpc wide and with a 
mass resolution of nip = 1.6 x 10^h~^ M©, the closest 
eight particles were iteratively collapsed, reducing the 
particle number to 128^ particles. These low resolution 
initial conditions were then evolved until z = 0, at which 
point eight clusters were selected in the mass range 1- 
3xl0^^/i^^ M0. All particles within two times the virial 
radius were then tracked back to their initial positions 
at z = 45, where they were regenerated to their origi- 
nal mass resolution and positions. These high resolution 
pockets are surrounded by a "buffer" zone with eight 
times the original mass resolution, which itself is nested 
in particles that are 64 times more massive than the par- 
ticles at the center of the cluster. These initial conditions 
were then re-simulated to z = 0, recording 63 outputs 
from z = 1.5 to z = so that At « 0.17 Gyrs. A sum- 
mary of the eight host halos is presented in Table [U and 
on quick inspection it should be immediately apparent 
that the eight hosts have widely varying masses and as- 
sembly histories. We calculate these quantities as follows: 
the virial radius is defined as the distance at which the 
average halo density drops below Phaioirvir) = AmrPb, 
where A.yir = 340 and ph is the cosmological background 
density. The virial mass is defined to be the mass in- 
side this radius. We calculate each host's age as the time 



elapsed since their fo rmation, which is defined, following 
iLacev fc Cold ()1993f ) , at the redshift where the halo first 
contains half of its present-day mass. 



Halo 


-Rvir 




■^form 


age 


A^sat(< fvir) 


# 1 


1.34 


2.82 


1.18 


8.37 


166 


# 2 


0.97 


1.05 


0.87 


7.17 


49 


# 3 


1.06 


1.38 


0.84 


7.01 


98 


# 4 


1.06 


1.38 


0.75 


6.57 


71 


# 5 


1.34 


2.80 


0.59 


5.65 


168 


# 6 


1.06 


1.39 


0.50 


5.06 


97 


# 7 


1.00 


1.16 


0.43 


4.52 


54 


# 8 


1.37 


3.00 


0.30 


3.42 


152 



TABLE 1 

Summary of the eight host dark matter halos at z = 0. 
Distance is measured in Mpc, mass in lO-'-^'/i"-'- Mq, and 
AGE IN Gyrs. Only satellites with more than 200 particles 

ARE tallied in THE LAST COLUMN. 



2.2. Identifying Substructure 

Simulation outputs merely tell us what the particle 
spatial and kinetic distribution is at each redshift. They 
give us no information about particle assignment or halo 
identity - which particle belongs to which halo? There is 
no unique answer to this question, mainly because there 
are many different ways in which a halo can be defined. 

A number of sophisticated algorithms ha ve been de- 
velop ed t o locate hal os within sim ulations (iDavis et al. 

19851: iFrenk et al.] [l988l : jBertschin ger fc Gelt 

19911 ISuto et all [l99l rv\^ber g eta ] [19971 : 
Klvpin fc HoltzmanI Il997f ). They face many chal- 
lenges: the dynamic environment of cosmological 
simulations blurs halo boundaries, and halos are con- 
tinually undergoing mergers or being stripped within a 
host potential, making it impossible to clearly define 
a halo edge. Furthermore, most of these do a poor 
job at finding substructure in very dense background 
regions, and although nearly all algorithms now use 
kinetic information to remove gravitationally unbound 
particles, they are generally not too concerned with 
background contamination, which can be safely ignored 
for most applications. 

Unfortunately, these issues are especially problematic 
for our analysis: once the substructure halos cross the 
virial radius of the host cluster, contrast is lost, and the 
particle background from the host becomes very signifi- 
cant. If we were to mistakenly assign background parti- 
cles from the cluster to our substructure halos, this could 
mimic the radial alignment effect we are looking for, since 
cluster particles themselves are radially distributed. We 
solve this problem by finding the substructure halos early 
on, at the formation redshift {zform) for each host. At 
these early times, the hosts are still starting to assemble 
and halos are not as clustered, and therefore much easier 
to identify. Once the halos have been found, their indi- 
vidual particle distributions can then be tracked forward 
in time through any environment, even into the dens- 
est cores of clusters, without suffering from background 
contamination. 

A more detailed description of our halo finding and 
tracking methods can be found in iGill et al.l ()2004[ ). so 
we provide here only a brief summary. We find and trun- 
cate all the halos in our simulation volume at Zform 
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using the AMIGA Halo Finder (AH F), the successo r 
of the MLAPM Halo Finder (MHF) (|Gill et al.l[200l . 
AHF uses the adaptive grids of AMIGA to locate ha- 
los within the simulation. AMIGA'S adaptive refinement 
meshes follow the density distribution by construction. 
The grid structure naturally "surrounds" the halos, as 
the halos are simply manifestations of overdensities. As 
AMIGA'S grids are adaptive it constructs a series of em- 
bedded grids, the higher refinement grids being subsets 
of grids on lower refinement levels. AHF takes this hi- 
erarchy of nested isolated grids and constructs a "grid 
tree". Within that tree, each branch represents a halo, 
thus identifying halos, sub-halos, sub-sub-halos and so 
on. 

Once we have found all the halos and sub-halos in our 
simulations at this redshift, we can start tracking their 
particle distributions through time. The main disadvan- 
tage of this method is that any subsequent accretion (af- 
ter Zform) onto the halos will, by design, be ignored. 
This seems a reasonable compromise - halos within halos 
travel through their environment too quickly to accrete 
significant amounts of particles and we assume that any 
particles acquired before the halo enters the host settle 
into the potential well isotropically, so that we are still 
obtaining a fair sample of the shape of the halos by only 
including the particles that were present at Zform- 

At each time step we look at the distribution of parti- 
cles for each halo and, after recalculating their center of 
mass, we check if each particle is still bound to the halo. 
This is an iterative process: starting at the center of the 
halo and moving outwards, we calculate each particle's 
kinetic and potential energy in the halo's reference frame 
and remove all particles that have velocities, v > bv^sc, 
where b — 1.5 is the bound factor, and the only free pa- 
rameter in our algorithm. We repeat the process until 
no further particles are removed or a minimum number 
{Np = 200, c.f. §2.3) of particles has been reached. Parti- 
cles that are determined to be unbound are subsequently 
ignored. This is a completely effective way of removing 
the cluster background, as particles that do not belong to 
the substructure halo will be quickly left behind. It also 
allows us to track debris being stripped off the subhalos 
as they orbit inside the cluster. 

When all unbound particles have been removed, we 
fit an NFW distribution to the radial profile of the re- 
maining particles. We define the halo's radius as the 
distance at which the average halo density drops below 
Phaioirvir) = ^vir{z)pb{z), whcrc Ay„{z) is the virial 
overdensity at that redshift, and discard any particles 
that lie outside this limit. However, this radius is almost 
never reached in the case of substructure, in which case 
the radius of the halo is defined as the distance to the 
furthest bound particle. Once we have determined which 
particles belong to which halo at each timestep, we are 
ready to measure their shapes. 

2.3. Shape Measurements 

How can we condense a three dimensional particle dis- 
tribution into a few simple parameters describing its 
shape? With no prior knowledge of how the particles 
are distributed this is a difficult task. However, ha- 
los produced in dark matter cosmologic al simulations 
seem to follow a universal density profile ([Navarro et al.l 
Il996l) . and are generally well fit by triaxial ellipsoids 



(|Frenk et al.lll988t lAllgood et al.ll2006f) . The simplest 
way to do this is to calculate the inertia tensor of the 
distribution, Ijk = '^i'>TT'ifi,jri^k, which is then diago- 
nalized to find the principal axes of the halo. However, 
this procedure is not ideal, since it weights particles by 

, and therefore results in a shape measurement that is 
overly biased by the outlying parti cle distribution. 

A better measure ()Gerhardlll983[ ) is the reduced inertia 
tensor: 

Iok = l^m,^^. (1) 

i » 

which weights particles equally regardless of their dis- 
tance to the center of the halo, using only the directional 
information of halo particles to calculate their shapes. 
The eigenvectors and eigenvalues of this reduced form 
of the inertia tensor give us the principal axes of the 
halo and a measure of their relative lengths {b/a,c/a), 
although the latter are substantially overestimated, as 
we shall see. 

The main source of uncertainty in determining the 
shapes of our halos is the small number of particles we 
sample their potentials with. We want to characterize 
halo alignments over as wide a mass range as possible, 
so we want to know what the minimum number of par- 
ticles is that will still give us a reliable measure of a 
halo's shape. The ability to determine the orientation of 
a halo's major axis also depends strongly on the value of 
b/a: An oblate halo with b k a, will be almost degenerate 
in its major/intermediate axis orientations. 

In order to address these questions, we generated a set 
of fake triaxial NFW halos with varying numbers of par- 
ticles {Np) and intermediate-to-major axis ratios {b/a) 
and fed them through our pipeline. For each value of 
Np and b/a we performed 100 random realizations of an 
NFW halo and calculated the angle between the major 
axis direction measured and that which was input. The 
dispersion in these values, 6acc, is then a good estimate of 
the accuracy of our measurement. The minor-to-major 
axis ratio (c/a) does not appear to affect the determina- 
tion of the major axis direction, and the results presented 
in figure [1] are therefore only for prolate halos with b — c. 

As expected, our accuracy depends very strongly on 
the number of particles sampled - the points on the left 
panel are well fit by a relation of the form: 6acc oc N"'^-^'^. 
We want a compromise between individual halo accuracy 
and sample size - at values of Np < 200, Oacc increases 
rapidly, and we pick this, somewhat arbitrarily, for our 
lower limit on Np. If b/a = 0.8, our measurements of 
these halos would be accurate by 6acc ~ 10°. However, 
9acc also depends strongly on b/a. When Np = 200, 
b/a < 0.8 is required to maintain the 10° error, with an 
increase in &/a leading to a rapid decrease in accuracy. 
Figure [1] refers to the input values of b/a. In fact, the 
measured ellipticities are much higher, although the two 
are tightly correlated: {b/a)^nput = {b/a)°Jelsured- We 
place an upper limit of 0.8 on the intrinsic axis ratios, 
which translates to a limit on the measured values of 
b/a < 0.9. 

With our limits in place for the minimum number of 
particles and maximum axis ratios, we are ready to start 
analyzing our results. It is worth noting, however, that 
both these error sources would bias our shapes randomly: 
there is no preferred direction that will be selected if the 
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halos are under-sampled or too spherical. This in turn 
implies that the results on alignment presented in the 
next section are, if anything, conservative. 

3. RESULTS 
3.1. Alignment at z — 

The quantity we will focus on is the angle, 0, between 
the major axis of each halo and the vector connecting 
the halo to the center of the host. If halos are oriented 
randomly in space, the cosine of (j) will be uniformly dis- 
tributed between and 1, with a mean value, (coS(/>), 
of 0.5. When cos0 ~ 1 the halo is pointing toward the 
host center, whereas when cos « it is aligned tangen- 
tially to it, so that a value of (cos0) > 0.5 implies an 
overall tendency for radial alignment. The standard er- 
ror on (cos0) is (T(cos0) = c/v^, where N is the sample 
size and a is its standard deviation. We show in Figure 
[2] a histogram of cos for all halos within 2rmr of each 
of the eight hosts. It is immediately apparent that our 
distribution is inconsistent with isotropy at a very high 
significance level: (cos0) = 0.66 ± 0.01, with most halos 
pointing toward the center of the host halo. 

While it is clear that the results of figure [2] confirm pre- 
vious observational reports of radial alignment, a precise 
quantitative comparison is rather difficult, and we de- 
fer this discussion to Nevertheless, much can be 
learned from a qualitative study of the effect's behavior 
and correlation with individual (and host) halo proper- 
ties. Figure [3] shows the same histogram as in Figure [2] 
but now for two separate halo populations, segregated 
by mass. There does not appear to be a significant dis- 
tinction between the two populations. This tells us not 
only that the alignment effect is mass independent, but 
also confirms the experiments in §2.3 that show that res- 
olution effects are unimportant in the lowest mass halos 
considered in our analysis {Np > 200). 

We can also study how the effect depends on extrin- 
sic characteristics of the halo, e.g. the distance to the 
center of the host, or the host mass. We searched for 
correlations with different global host properties such as 
mass and age, and found none. The alignment mecha- 
nism appears to be universal, in that it is present with 
approximately the same strength in hosts with widely 
varying mass, formation times and assembly histories. 
Thi s surprising resul t is als o seen in observational stud- 
ies: iPereira fc KuhnI ()2005[ ) found no correlation of the 
alignment strength with the dynamical state of the clus- 
ters inferred from their x-ray morphologies. 

3.2. Dependence on Distance to Cluster Center 

Figure [3] shows the dependence of the effect on the 
distance from the cluster center. All halos at redshifts 
z < Zforrn SLie iucludcd in this analysis, in order to 
enhance the overall signal. The behavior appears very 
smooth: (cos (f>) rises gradually as the host is approached, 
peaks slightly past its virial radius, and then decreases 
again toward the center. 

It is striking that already at a distance of three virial 
radii there is a small, consistent, tendency for radial 
alignment. At this distance, how can the halo already 
"feel" the presence of the host? This is easily under- 
stood once we consider that clusters form at the inter- 
section of filaments, and hence that most filaments will 



be radially aligned with respect to their nearest clusters. 
If there is a primordial alignment of halos with respect 
to the filaments in which they form, then even at large 
distances this will appear as a radial alignment in our 
analysis. T his type of primordial alignment at large radii 
was seen bv lAragon-Calvo et al.l (2Q07) in their study of 
filamentary structures, where they found similar values 
for (cos0) (their figure 2e). 

The main focus of this paper, however, is what hap- 
pens closer to the host. As the halo falls in, the am- 
plitude of the alignment increases dramatically, reach- 
ing a peak of (coscj)) — 0.72 at about one-half of the 
virial radius of the cluster, before decreasing again grad- 
ually inside the core. The increase of the alignment 
with decreasing dis tance matches the behaviour found by 
iFaltenbacher et al.l ( |2007 ) in their study of SDSS groups, 
and is to be expected if the effect is caused by the tidal 
field of the host, but the dip at small radii, r < 0.3r„ir, 
has not yet been observed. This is most likely due to 
the severe projection effects that dominate the cores of 
observed clusters. 

What causes this behavior? It appears that the align- 
ment that is set-up in the infall regions is being disrupted 
in the inner regions of the cluster. What causes the dis- 
ruption? Is this primarily a spatial effect caused by the 
environment of the cluster core, or a temporal one, given 
that galaxies closest to the center have been in the cluster 
environment for longer? And what produced the align- 
ment in the first place? The best way to answer these 
questions is to take advantage of the extra dimension pro- 
vided by simulations, and explore the evolution of this 
effect with time. 

3.3. Evolution with Redshift 

The evolution of the alignment with redshift is plot- 
ted in figure [5] for each of the eight clusters indepen- 
dently. Perhaps the most striking feature of this plot 
is the self-similarity of the different curves. Every clus- 
ter appears to go through exactly the same evolution, 
regardless of size or formation time, such that at 2; = 
they are practically indistinguishable, as described in the 
previous section. Figure [5] also reveals that the clusters 
evolve monotonically, with the strength of the effect in- 
creasing steadily since the formation time of each cluster 
to the present day. 

Whatever the source of the disruption at the cluster 
cores, it is seemingly not strong enough to dilute the over- 
all alignment signal. There are two possible explanations: 
It could be that, even though alignment is disrupted once 
the halo reaches the core of the host, the constant infall 
of pristinely aligned halos results in an overall increase of 
the average alignment per host. Alternatively, the mis- 
alignment seen at the cores could be short-lived - a fea- 
ture of each halo's orbital motion through the potential 
of the host. 

Distinguishing between these two alternatives requires 
a different approach: we need to track halos on their way 
toward the cluster, and then trace their orbits inside the 
virial radius of the host. 

3.4. Evolution with Orbital Phase 

Figure [H] shows the alignment evolution stacked for all 
halos throughout their orbits. Initially halos are tracked 
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relative to the amount of time (in Gyrs) remaining un- 
til they cross the virial radius of the host for the first 
time. Once they cross this threshold, halo orbital times 
are normalized at each passage through pericenter and 
apo center. 

We again detect a small alignment at large distances 
from the cluster , which we believe is evidence for a pri- 
mordial alignment along filaments as discussed in ^3.21 
As the host is approached the signal increases signifi- 
cantly, peaking just before the first pericentric passage, 
and then a periodic oscillation ensues, which follows the 
halo's orbital period closely. On average, the tendency 
for alignment is much larger within the host than before, 
although the alignment tendency changes dramatically 
with orbital phase, ft follows that the dip observed near 
the cluster cores in figure 3] is in fact a result of the mis- 
alignment observed at pericenter, and, most importantly, 
that it is not disruptive, since the alignment tendency is 
restored well before the next apocenter is reached. In 
fact, the alignment is quite constant throughout the rest 
of the orbit and seems to increase slightly at each pas- 
sage. This evidence points to a stable dynamical effect 
that is set-up as the halo orbits around the cluster. 

Further insights can be obtained by exploring the ori- 
entations of the halos with respect to their orbits. We 
define a new angle, f3, as the angle between each halo's 
major axis and the halo's velocity, and plot the mean 
value of its cosine for all halos vs. orbital phase in Figure 
[71 The similarities in behavior between radial and orbital 
alignment at large distances are simply a consequence of 
the radial nature of the orbits themselves - halos form 
and travel along filaments toward the intersecting nodes 
where clusters reside. In fact, even inside the hosts, or- 
bits are quite eccentric, with an average apocentric to 
pericentric distance ratio of 4 : 1. 

Could it be, then, that the radial alignment we observe 
within the virial radius is just a tendency for halos to 
be aligned along their orbits coupled with the fact that 
orbits are, on average, quite radial? Figure [7] tells us 
that this is not the case: once inside the cluster, we find 
that the orbital alignment is also correlated with orbital 
phase, but whereas the radial alignment is almost in- 
stantly recovered after pericenter, the orbital alignment 
increases again much more slowly, and only after reach- 
ing the next apocenter. This asymmetry around peri- 
center seems, at first, surprising, but, as will be shown 
in the following section, follows as a natural consequence 
of tidal torquing by the cluster potential throughout the 
halo's orbit. 

4. DISCUSSION 

4.1. Tidal Torquing as a Mechanism for Alignment 

Once radial and orbital alignment information is com- 
bined, a clearer picture emerges of what is going on in- 
side these clusters. As the halo approaches pericenter 
along an eccentric orbit, it is continually torqucd along 
the direction of the potential gradient, i.e. halos tend to 
point toward the host center, and, because their orbits 
are fairly eccentric, also along their orbital direction. At 
pericenter, the halo is moving too fast for the torquing 
to be completely effective, which causes the dip in radial 
alignment. It is nevertheless enough to torque the halo 
away from its orbital direction and back toward the clus- 



ter center, in a figure rotation that is co-planar with its 
orbital rotation and in the same direction. The radial 
alignment is quickly reinstated, but orbital alignment is 
lost as the halo progresses towards apocenter. Steady 
torquing throughout the orbit keeps halos oriented to- 
ward the cluster center and away from the direction of 
their orbits until after the apocentric passage, where or- 
bital alignment increases steadily towards pericenter, and 
a new cycle begins. Figure [8] illustrates this behaviour 
with a sketch of a halo's rotation as it orbits around the 
cluster. 

If halo orbits were circular, halos would quickly become 
tidally locked and maintain radial alignment throughout 
their orbits. In reality, their orbits are quite eccentric, 
and their orbital speed varies significantly. Halos do not 
react to the tidal torquing quickly enough through the 
pericentric passage, and the narrow dips observed are 
the result. In fact, idealized numerical experiments in- 
volving a single halo in a circular orbit around a static 
host invariably lead to tidal locking of the halo, although 
the time required for locking varies significantly with the 
original orientation of the halo (C. M. Simpson & K. V. 
Johnston, private communication). Interestingly, for ha- 
los that start out already pointing toward the host cen- 
ter, the time required is rather short, of the order of an 
orbital period or less. 

Further support for this tidal torquing hypothesis is 
shown in figure O Although we believe our shape mea- 
surements to be robust to random outliers, it is possible 
that strongly distorted outer shells, caused, e.g., by tidal 
stripping, could significantly bias the result. As a test, 
we apply four different particle cuts to each of our halos 
by varying the boundedness criteria on the particle ve- 
locities: instead of throwing out all particles for which 
V > bvesc, where b = 1.5, we exclude alternately particles 
that have velocities greater than 1, 0.75 and 0.5 times the 
escape velocity. For the most conservative criteria, which 
only retains particles that have velocities v < 0.5vesc, 
more than 70% of the particles are discarded, and we 
are only probing the very bound cores of the halos. Fig- 
ure [9] makes clear that stripping cannot possibly be the 
sole cause of the alignment effect, since even the most 
conservative cut shows significant alignment. 

Nevertheless, a trend is observed, in that the most 
bound particles show slightly less tendency for alignment 
overall. This could be the result of tidal stripping in the 
outer layers, but more likely it is a simple statistical ef- 
fect: Because we only consider halos with Np > 200 in 
this analysis, as we progressively exclude more particles 
from the halo with decreasing b, some halos fall below 
this limit and are consequently ignored. Hence a decrease 
in b implies a smaller number of halos in each sampled 
bin, which reduces the signal-to-noise, and brings (cos (f)) 
closer to 0.5. Despite this trend, the conclusion remains 
that halo shapes are not significantly warped by tidal 
stripping, and that tidal torquing of the entire halo is a 
better explanation for the effect. 

A number of early numerical and analytical studies 
support the importanc e of tidal torques within clusters . 
iMiller fc SmithI ()1982f ) performed a set of numerical ex- 
periments on a rotating bar in an external force field and 
observed tidal braking of the rotation, with a rate that 
was inversely related to the square of the cluster crossing 
time. More recently, numerical A^-body experiments by 
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ICiotti &: Duttal (|1994[ ) showed that the time required for 
the ahgnment of a prolate galaxy with the tidal field of 
a cluster is much shorter than the Hubble time, and on 
the order of a few times the galaxy' s intrinsic dynami- 
cal tim e. Using a different approach, Usami fc Fujimotol 
(|1997t l studied tidal effects on gaseous ellipsoids orbiting 
in a central potential analytically, predicting that galax- 
ies in eccentric orbits should have their long-axis trapped 
toward the direction of the radius vector of the cluster. 

While this paper was being written, two studies were 
publish ed on halo aU g nmen ts that describe similar re- 
sults. iKuhlen et al.l (|2007D studied the alignment of 
substructure around a Milky Way type halo using the 
Via Lactea simulation, and observed a radial alignment 
tend ency that is preserved th roughout the halos' or- 
bits. iFaltenbacher et al] (|2007[ ) looked at several differ- 
ent types of alignment in a set of dark matter hosts at 
z = 0, finding similar levels of radial alignment that in- 
crease with decreasing distance to the host. 

4.2. A Comparison with Observations 

The results of ^3] certainly seem to substantiate the ob- 
servational evidence for radial alignment of cluster galax- 
ies. A quantitative comparison, however, is not easily 
made. In order to properly "observe" these simulations, 
semi-analytic models of galaxy formation are required to 
extrapolate from the dark matter halos to the luminous 
components embedded within. These then need to be 
projected, interlopers and survey limits accounted for, 
and the resulting image fed through traditional source 
extraction and isophotal analysis pipelines. This is a 
laborious procedure and it cannot yet provide accurate 
results, since galaxy formation models are still largely 
unconstrained in a crucial parameter: the alignment be- 
tween luminous and dark matter. 

Observationally, studies of the alignment and relative 
ellipticity of the two components are currently only pos- 
sible for gravitational l e ns ga laxies, a very rare class of 
objects. iKeeton et al.l ()1998[ ) analysed a sample of 17 
lenses, mostly isolated early-types, and found that the 
luminous component of the lens generally aligns with its 
inner halo to < 10°. In order to probe the shapes of the 
halos to larger radii, stacked galaxy-galaxy weak lensing 
studies are needed. These are just now becoming feasible, 
and preliminary results appear somewhat contradictory 
iHockstra ct al. 2004; Mandclbaum et al. 2006b). 

Most theoretical studies have concentrated on the for- 
mation of disk galaxies and their angular momentum, 
where some misalignment betwee n baryonic and dark 
matter spin is commonly seen (e.g. 'van den Bosc h et al.l 
{2002)). On the other hand, Baihn ct al. (2005) find that 
the orientations of simulated halos and their embedded 
disks are largely uncorrelated at large radii, and almost 
perfectly ahgned at small (r < O.lr^i^). 

The current uncertainty in this parameter makes it 
impossible to accurately predict the orientation of the 
galaxies that would populate our halos. However, the 
results presented in this paper suggest a gravitational 
origin for the alignment mechanism, and it is therefore 
reasonable to expect that the two components should 
react similarly to it. Furthermore, the tidal torquing 
within clusters is so effective that the halos appear to 
"forget" their original orientations before a single orbit 
is completed, which renders the original alignment be- 



tween light and dark matter relatively unimportant. 

We therefore compare the dark matt er alignment di- 
rectly with the galaxy observations of iPereira fc Kuhnl 
(|2005f ). We project each halo's dark matter particles 
along the three spatial dimensions in our simulation and 
compute the 2D inertia tensor of their projected distribu- 
tion. The angle between the halo's 2D major axis and its 
projected separation from the cluster center can then be 
measured. We include in this sample all galaxies within 
2 virial radii of the cluster center - interlopers are not 
accounted for, since the SDSS galaxies we are comparing 
our results to are all spectroscopically confirmed cluster 
members. Figure [TU] shows the results of this 2D analy- 
sis and compares them with the SDSS observations. We 
plot all three independent projections and note that the 
dispersion in their values should give us a fair estimate of 
the error introduced by the projection procedure itself. 

The dark matter alignment is much stronger than 
that observed: (6I)^„,„ = 34°. 5 ± 0°.9 whereas (0)^^, = 
42°. 79 ± 0°.55. This must be in some part a reflection of 
how much harder it is to measure accurate galaxy posi- 
tion angles on an survey image than for a well-resolved 
halo in a cosmological simulation, where dynamical in- 
formation allows for a much cleaner background removal. 
Nonetheless, the dilution caused by this measurement 
noise cannot wholly account for the significant difference 
in radial alignment between the two components. Given 
the nature of the alignment mechanism established in 
it is perhaps not too surprising that dark matter ha- 
los are more strongly aligned. One would naively expect 
the dark matter halos to be more easily torqued, given 
that they are much more extended (providing a longer 
lever) and have generally lower spins ( and therefore less 
gyroscopic resistance) than their luminous counterparts. 

4.3. Possible Consequences of Tidal Torquing in 
Clusters 

Halo alignments have traditionally been studied ei- 
ther as a probe of their formation history, or as a con- 
taminant to weak lensing studies. Now that we have 
established that the leading mechanism behind halo 
alignments within clusters is a dynamical effect present 
throughout their lifetime, it is interesting to speculate on 
what possible evolutionary consequences this mechanism 
might have for the halos affected. 

Figure [5] shows us that at each point in the orbit, the 
torque acts to rotate the halo away from its orbital di- 
rection, which necessarily results in a deceleration of the 
halo's orbital motion, inducing orbital decay. The ha- 
los analysed in this stu dy do indeed sho w a tendency for 
orbital circularization: I Gill et al] ()2004[ ) showed that ha- 
los with more pericentric passages have smaller orbital 
eccentricites. They argued that dynamical friction was 
not a likely cause, and tentatively ascribed the effect to 
the growth of the host halo instead. While the velocity 
dispersion of the satellites is seen to depend on the host 
halo mass, it seems possible that at least part of the or- 
bital decay observed is a natural result of the constant 
torquing throughout the halo's orbit. This is an interest- 
ing prospect, since an extra source of orbital decay could 
potentially help solve the outstanding problem of cD for- 
mation in massive clusters, as well as alleviate some un- 
resolved discrepancies between observed satellite popula- 
tions and the generally low efficacy of dynamical friction 
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predic t ed by numeri cal studies (e.g. iHashimoto et alj 
(|2003l ): lTaffoni et all (|2003l )). We are currently investi- 
gating the importance of this induced orbital decay, and 
this will be the subject of a future paper. 

Another possible consequence of the strong torquing of 
dark matter halos within hosts is the possibility of disk 
warping. Because of their high angular momentum, disks 
will naturally resist tidal torquing more effectively than 
the surrounding dark matter halo, which will introduce 
a misalignment between the halo and the disk. Even 
though recent studies of (isolated) disk-halo alignments 
show that t heir orientations a re largely uncorrelated at 
large radii (jBailin et al.ll2005l ). the same is not true for 
the inner halos (r < 0.1r„ir), where the rotational axis of 
the disk is seen to lie very close to the minor axis of the 
inner halo. We have shown that tidal torquing affects all 
particles in the halos, even the most bound, so it is not 
unreasonable to expect that the inner shells should also 
feel these torques. The question then remains whether 
the disk within will align itself accordingly, or whether 
the misalignment could be a possible cause of warping of 
the disk, but this will also require further study. 

5. CONCLUSIONS 

There is growing observational evidence that a satel- 
lite's major axis is preferentially aligned with the ra- 
dial vector linking the satellite to its host. This ten- 
dency for satellites to point at their hosts has been seen 
on both cluster a nd group scales (Pereira & Kuhn 200|; 
lAgustsson fc BrainerdI 12006) . Motivated by this result, 
we have used a suite of cosmological N-body simulations 
to investigate the alignment between satellite and host 
dark matter halos. 

We take particular care to separate satellite and cluster 
particles using a combined halo find er plus tra,cker. In 
this method, an adaptive halo finder ()Gill et al.ll2004D is 
used to initially identify a set of satellite sub-halos that 
we subsequently track as they enter and orbit the clus- 
ter, removing particles as they become unbound. The 
advantage of this approach is that we can be sure to 
use only genuine sub-halo particles and exclude "back- 
ground" cluster particles that might bias our shape mea- 
surements. We then use the reduced inertia tensor to 
measure the shapes and orientations of all sub-halos 
which end up inside the virial radii of a set of eight 
simulated clusters. We highlight here the main results 
obtained from this analysis: 

• The satellites in the simulations show a strong ten- 
dency to point toward the cluster center. The 
mean cosine of the angle between the major axis 
of each halo and the cluster center is (cost/)) = 
0.66 ± 0.01, where an isotropic distribution would 
have (cos(/)) — 0.5. This tendency for alignment 
is found for all clusters at all redshifts analyzed, 
and does not appear to depend on the mass of the 
cluster or the satellite. 

• The amplitude of the alignment is a strong func- 



tion of radius, with a small but significant effect 
extending out to many virial radii from the clus- 
ter. This signal, which has bee n seen in previous 
work ()Arag6n-Calvo et al.ll2007[ ). is most likely left 
over from the primordial imprint of the surrounding 
large-scale structure and can be ascribed to tidal 
torques exerted at early times, when the cluster- 
size per turbations were just turning around (e.g., 
lPeeblesl [r969'). 

• Closer to the cluster center, within 1-2 virial radii, 
the amplitude of the alignment increases dramati- 
cally to a peak of (cost/)) — 0.72 at about one-half 
of the virial radius, and then falls slowly closer to 
the cluster center. 

• When examined as a function of orbital phase for a 
given satellite, we find that the alignment increases 
rapidly as the satellite falls into the cluster for the 
first time and remains high after that, except for 
a short period during pericenter passage, when it 
dips precipitously. It is this short-lived dip which 
gives rise to the decrease in (cost/)) close to the 
cluster center. 

Based on these results, we conclude that the strong 
alignment seen at small radius — within two virial radii 
— is due to tidal torquing by the cluster halo. The idea 
is very simple - the galaxy is only in a stable equilibrium 
if it is pointing at the cluster center; otherwise there is a 
net torque which acts to rotate the galaxy towards this 
equilibrium point. We demonstrate that the alignment is 
seen both in the outer and inner parts of the satellite, in- 
dicating that it is not due to some process (such as tidal 
stripping) which impacts only the outer, poorly bound, 
part of the sub-halo. We also briefly review previous lit- 
erature which has investigated the impact of tidal torques 
on coUisionless systems using analytic approximations or 
idealized simulations, and find that the expected ampli- 
tude and timescale is sufficient to produce the alignments 
we see. 

Although we study only dark-matter halos, we expect 
this effect to extend to the luminous part of galaxies, as 
observations seem to indicate. This will have an obser- 
vational impact on weak lensing studies and may also 
modify the distribution of stars in a satellite, as well 
as the satellite's orbital properties. We will investigate 
these possibilities in future work. 

The simulations presented in this paper were carried 
out on the Beowulf cluster at the Centre for Astrophysics 
& Supercomputing, Swinburne University. We would like 
to thank Jeff Kuhn, Kathryn Johnston and Christine 
Simpson for helpful discussions. Greg Bryan acknowl- 
edges support from NSF grants AST-05-07161, AST-05- 
47823, and AST-06-06959, as weU as the National Center 
for Supercomputing Applications. 
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Fig. 1. Results for "fake" NFW halo analysis. Left: Error in major axis orientation vs. number of particles in halo. The vertical 
dashed Hue represents the lower limit on N for the lialos analysed in this paper, whereas the horizontal lines show the minimum accuracy 
that will be tolerated. Right: The error in major axis orientation vs. the intermediate-to-major axis ratio. Errors are shown for halos 
with 200 particles (as dots) and with 1000 particles (as crosses). 
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Fig. 6. — Average radial alignment vs. orbital phase for all sub- halos. The horizontal axis represents time - initially in Gyrs before 
crossing host R„ir a-nd then in subsequent Pericentric (P) and Apocentric (A) passages. 
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Fig. 7. — Average orbital alignment vs. orbital phase for all sub- halos. The horizontal axis represents time - initially in Gyrs before 
crossing host R„ir and then in subsequent Pericentric (P) and Apocentric (A) passages. 




Fig. 8. — A sketch of a representative halo orbit around its host. The center of mass of the host is represented by a star. The radial 
direction is drawn as a red vector, cyan indicates the orbital direction, and the purple vector is the direction of the major axis of the halo. 
Purple and red vectors are generally very close, apart from a short mismatch at pericenter caused by the high orbital velocities. Red and 
cyan vectors are close before pericenter, but almost orthogonal to each other in the second part of the orbit. 



13 




0.25 



0.20 



0.15 



o 

X 

o 

c 

o 



g 0.10 




0.05 - 



0.00 

40 60 
( degrees ) 

Fig. 10. — 2D radial angle distribution at 2 = for all hales within 2ryir of the cluster centers, projected along the three independent 
spatial dimensions (filled circles). Also plotted for comparison is the observed angle distribution from iPercira fc Kulin [2005.) . with 
Poissonian error bars. An isotropic distribution would follow the horizontal dashed line, with (6) = 45° 



